Rotationally driven gas of inelastic rough spheres 
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We study a two-dimensional gas of inelastic rough spheres, driven on the rotational degrees of 
freedom. Numerical simulations are compared to mean-field (MF) predictions with surprisingly 
good agreement for strong coupling of rotational and translational degrees of freedom - even for very 
strong dissipation in the translational degrees. Although the system is spatially homogeneous, the 
rotational velocity distribution is essentially Maxwellian. Surprisingly, the distribution of tangential 
velocities is strongly deviating from a Maxwellian. An interpretation of these results is proposed, 
as well as a setup for an experiment. 
PACS: 45.70, 47.50+d, 51.10.+y, 47.11.+j 
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Systems of hard spheres have a long-standing history 
as a basic model for gases, liquids, and e.g. glasses. When 
dissipation is added, one has the minimal model for gran- 
ular materials and it is only a small step to include also 
the rotational degrees of freedom via a tangential inter- 
action at contact. Granular materials belong to the fas- 
cinating world of non-linear, dissipative, non-equilibrium 
systems , whose interest is due to their practical im- 
portance and due to the theoretical challenges they rep- 
resent. Granular media are collections of macroscopic 
particles with rough surfaces and dissipative, frictional 
interactions. Numerical simulations are an established 
tool to complement advanced theoretical approaches and 
difficult experimental studies. 

In order to study systems of rough spheres, kinetic 
theories have been extended to (weak) dissipation and 
friction [^|J^]. Alternative, more recent approaches are 
based on a pscudo-Liouville operator formalism |7] [R| 
and are less general in the sense that they assume homo- 
geneity and Maxwellian velocity distributions in order to 
arrive at a mean field (MF) description of systems with 
rotational degrees of freedom. Either the system is left 
undisturbed [[7[-|| and thus cools continuously, or a "driv- 
ing" can be applied, i.e. energy is fed into the system in 
order to reach a steady-state situation. 

The typical driving of a granular material, in both 
experiment and simulation, can be realized by moving 
walls |Q] which lead to rather localized input of energy. 
Alternatively, the system can be driven by a global ho- 
mogeneous, random energy source in different variations 
P|-^5[. Depending on the experimental setup, energy 
can be given either to translational degrees of freedom or 
to the rotational ones, or to both. The first case catchcd 
most of the attention - reason enough to change the focus 
and feed rotational energy instead of translational. In the 
experiment, translational energy input was applied for 
special boundary conditions and a variety of interesting 
experimental results were obtained just recently [^G| p(j[ . 
One can obtain a gas and a liquid state, together with 
dense, solid-like clusters which form due to dissipation. 



The dynamics of the system is usually assumed to be 
dominated by two-particle collisions which are modeled 
by their asymptotic states: A collision is characterized 
by the velocities before and after the contact, and the 
contact is assumed to be instantaneous. In the simplest 
model, one describes inelastic collisions by a normal resti- 
tution coefficient r only, i.e. the negative ratio between 
the normal velocities after and before the collision. How- 
ever, since surface roughness and friction are important 
||7|- |To| , pl 22 1, one should allow for an exchange of trans- 
lational and rotational energy. In the simplest approach 
[p|,|7[|2l| , surface roughness is accounted for by a constant 
tangential restitution coefficient 7*t, which is defined in 
analogy to r in the tangential direction. A more realis- 
tic friction law involves the Coulomb friction coefficient 
i0j23- 



-25||, so that the tangential restitution will de- 
pend on the collision angle. Constant tangential restitu- 
tion is recovered in the limit of perfect friction. 

In this Letter, we will focus on a system of such per- 
fectly rough particles, where only the rotational degrees 
of freedom are coupled to a homogeneous driving. Such a 
situation could correspond, for example, to a gas of rough 
magnetic particles subject to a rapidly varying, homoge- 
neous, magnetic field. Besides a possible experimental 
application, we believe that this study is interesting in 
itself, since a correct modeling of the driving mechanism 
is of great importance for a theory of granular gases to 
describe realistic experimental situations. 

The model consists of N three-dimensional spheres 
with radius a and mass M , interacting via a hard-core po- 
tential and confined to a 2D plane of linear extension L, 
with periodic boundary conditions. The degrees of free- 
dom are the positions Ti(t) the translational velocities 
Vi(t) and the rotational velocities uJi(t) for each sphere 
numbered by i = 1, . . . , N. When two particles 1 and 2 
collide, their velocities after collision are related to the ve- 
locities before collision, through a collision matrix which 
is derived from the linear and angular momentum conser- 
vation laws, energy/dissipation balance. The magnitude 



of dissipation is proportional to the quantities 1 — r and 



1 



1 — rf , while the strength of the coupling between rota- 
tional and translational motion is connected to 1 + r t , 
where the normal restitution r varies between 1 (elastic) 
and (inelastic) and the tangential restitution r t varies 
between —1 (smooth) and +1 (rough), corresponding to 
zero and maximum coupling, respectively H|t[|2~^] . For a 
typical steady-state configuration (energy input is speci- 
fied below), see Fig. |l|. 
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FIG. 1. Snapshot of the particle distribution in the steady 
state for a system of N — 11025 particles, v = 0.34, r t = 1, 
and r n = 0.1. 

In order to feed energy, the system is agitated each 
time interval At = fJ-, with a driving rate /d r - Here, we 
will not apply driving frequencies much higher than the 
collision rate f2, but will use driving frequencies around 
100 s _1 , comparable to Q. This is rectified, since numer- 
ical checks with strongly different values of /dr lead to a 
similar behavior of the system even for driving frequen- 
cies lower than, but of the same order as Q, provided 
that a stationary state is reached. The translational ve- 
locity remains unchanged, but the angular velocity of 
particle i is modified at each time of agitation t so that 

a4(i) =0Ji(t) +riUjQ , 

where the prime on the left hand side indicates the value 
after the driving event. Due to the two-dimensionality 
of the system, we apply the driving force only to the z- 
direction, so that the scalar u is to be understood as the 
z-component of u). loq is a reference angular velocity (in 
this study we use ojo = 2.4 10 _4 s _1 ) which allows, with 
vq = auto = 2.4 10 _7 ms _1 , where a — 10~ 3 m, to define 
the dimcnsionless translational and rotational particle 
temperatures T tr = E tT /(NT ) and T rot = 2E mt / \NT ), 
with the translational energy Et T = (M/2) ^2i=i v ii the 
rotational energy E TOt = (qMa 2 /2) ^ i=1 uj 2 , and the ref- 
erence temperature To = Mv 2 . The variance of the un- 
corrected Gaussian random numbers (with zero mean) 
can now be interpreted as a dimensionless driving tem- 
perature Td r pT[ . The stochastic driving leads thus to 
an average rate of change of temperature 



AT mt /At = H dl , with H dl = / dr T ( 



dr 



(1) 



The starting point for our mean-field analysis is the 
theory of Huthmann and Zippclius M , for a freely cool- 



ing gas of infinitely rough particles, which was recently 
complemented by numerical simulations in 2D and 3D 
B and by studies of driven systems as well [jll). The 
main outcome of this approach is a set of coupled evo- 
lution equations for the translational and rotational MF 
temperatures T tr and T rot |?J which can be extended to 
describe arbitrary energy input (driving) |Tlj |. In the 
present study, given the random driving temperature T dr 
and an energy input rate /d r , as defined above, one just 
has to add the positive rate of change of rotational energy 
i?d r to the system of equations: 
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FIG. 2. Simulation (points) and theory (lines) for the pa- 
rameters v — 0.34, N = 11025, and r t — 1, plotted against r. 
(a) Equilibrium rotational temperature T ro t, normalized by 
the MF value T^ t (r = 0) at r = 0. (b) Ratio of equilibrium 
rotational and translational temperature R = T rot /T tr . 

with G — 8/ (^/itMa)vg2a(v), and the pair correlation 
function at contact g2a(v) = (1 — 7^/16)/(l — v) 2 in 



the approximation proposed by Henderson [ p6|j27| ] , de- 
pendent only on the volume fraction of the granular gas 
v = na 2 N /V . The constant coefficients in Eqs. (Q) and 
(|) are A = (1 - r 2 )/4 + 77(1 - T])/2, B = r 1 2 /(2q), and 
C = rj (1 — rj/q) /(2q), with the abbreviation r\ — rj(r t ) = 
q(l + r t )/(2q + 2), as derived in Ref. @, 

Setting to zero the temporal derivatives in Eqs. (^|) and 
(0) , one obtains the steady state properties of the system: 
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and 



T£ = TZ/R 



(4) 



with r = 2 {B/Ay' z (C - B 2 /A), and R = A/B. 
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In Fig. H we present equilibrium values of T rot , normal- 
ized by the MF value T r e 9 t (r = 0) at r = 0, and of the 
ratio R = T°^ t /T°^, as obtained from numerical simula- 
tions of a system of TV = 11025 particles, with volume 
fraction v = 0.34, r% = 1, and r ranging from 0.99 to 
10~ 4 . Surprisingly, the agreement with the MF predic- 
tion is very good, even for the lowest value r = 10~ 4 of 
the normal restitution, which corresponds to very strong 
dissipation, were the deviation from MF theory is of the 
order of only 10%. 
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FIG. 3. Steady state rotational (a) and translational (b) 
velocity distributions for TV = 11025, v = 0.34, r t = 1.0 and 
T = 0.1. A power law fit (dashed line) gives an exponent 
a — 1.92(6) for the rotational distribution and a = 1.41(6) 
for the translational distribution (see text for details). For 
comparison, a Maxwellian (full line) is plotted in (b). 

To give an example, if the system is driven on the 
translational degrees of freedom, the equilibrium temper- 
atures show deviations of 30 — 40% from MF predictions 
already for r = 0.6, sec [jfl). The snapshot in Fig. || 
shows the equilibrium particle distribution for r = 0.1 
and appears spatially homogeneous - apart from small 
density fluctuations not quantified here. 

In Fig. ^, we show the equilibrium rotational and 
translational velocity distributions for r = 0.1, with 
the other parameters as above. The rotational veloc- 
ity distribution is very near to a Maxwellian. A three 
parameter fit f(x) — Aexp(— B\(x — (x))/a\ a ) 1 where 

1/2 

o~ = \{(x — (x)) 2 )) , and x either equals u> or v, is 
plotted as dashed line in Fig. ||. The parameters (uj) 
and a are taken from the simulations, and the fit gives 
q = 1.92(6) for u>, while we obtain a = 1.41(6) for v. 



Rotational velocities are then characterized by good ho- 
mogeneization at low r, while the translational velocity 
distribution shows strong deviations from a Maxwellian. 
This deviation is due to the high dissipation. Numerical 
simulations with r = 0.99 give a Maxwellian distribution 
for both rotational and translational velocities. 

In order to check the role of the tangential restitu- 
tion, we show in Fig. || the equilibrium values of R with 
r = 0.1 and r< € [—1,1]. While for positive rt there is 
still good agreement with MF theory, strong deviations 
appear as r t — > — 1. Note that many realistic materials 
obey the relation r t ~ 0.4 [^4|, what renders our mean 
field approach still acceptable. 
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FIG. 4. Simulation (points) and theory (lines) results for 
R = Trot/Itr, with parameters v = 0.34, N = 11025, and 
r — 0.1, plotted against rt- 

The conclusions of our study are that the driving on 
the rotational degrees of freedom is able to keep the spa- 
tial homogeneity of the system up to very high dissipa- 
tion rates, for positive values of r t . This leads to a very 
good agreement of the equilibrium temperatures with the 
MF predictions. There are two possible reasons for this. 
From one side, the driving acts on rotations. Then, it 
cannot favorise collisions, since it does not increase the 
normal component of the relative velocity of the colliding 
particles. From the other side, the increase of rotational 
energy triggered by the driving leads to a shearing force 
between particles, which reduces density fluctuations and 
destroys velocity correlations. When r t — > —1, the agree- 
ment with MF is lost. To explain this result one has to 
remember that 1 + r t is a measure for the strength of the 
coupling. Not enough rotational energy is transferred to 
the translational degree, so that the randomization on 
collision does not take place. Thus, it is not surpris- 
ing that MF is no more valid in this very singular limit. 
Snapshots of the particle distribution for r = 0.1 and 
r t near to —1 (not displayed here) show indeed stronger 
density fluctuations in the system as reported in Fig. [j]. 

The translational velocity distribution exhibits strong 
deviations from MF prediction, in the homogeneous high 
dissipation regime, showing from one side that deviations 
from a Maxwellian are not necessarily related to cluster- 
ing, and from the other side that it is possible to have 
a good agreement of the second moment (the tempera- 
ture) of the velocity distribution with MF theory together 
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with a non Maxwellian velocity distribution. This poses 
a theoretical challenge, since recently proposed theories 
for horizontally driven granular gases ^8|, assume that 
clustering is responsible for fat tails in the velocity dis- 
tribution. 

Apparently the temperature of the system depends 
mainly on the energy balance relations, which depends 
indirectly on density fluctuations (density fluctuations in- 
fluence strongly both the frequency of collision and the 
rate of dissipated energy per collision) while higher mo- 
ments, and the overall shape of the velocity distributions 
are more sensible to other details. 

A difficult to realize, but probably good, experimental 
setup is the following. Each, extremely rough, granular 
sphere, contains a small (to reduce the effect of dipole- 
dipole interaction at collision) magnetic bar. The plane 
on which the spheres move should be extremely smooth, 
in order to avoid energy dissipation. Then, a spatially 
homogeneous magnetic pulses periodically spaced in time 
can be applied in the horizontal directions. This would 
be the magnetic analogon of the oscillating plane. If the 
magnetic field is really spatially homogeneous, the mag- 
netic dipoles of the spheres will receive angular momen- 
tum from the field, so only rotations are driven, and this 
angular momentum will be "randomized" by the colli- 
sions, if they are frequent enough, exactly like it happens 
for the kinetic energy coming from the oscillating plane 
usually employed in experiments. To reach an equilib- 
rium, it is necessary to give an initial translational veloc- 
ity to the particles. We are aware that such an experi- 
ment is extremely difficult to realize, but we hope there 
will be some experimentalist willing to try it. 

Summarizing, the main discovery of this work is that a 
dissipative gas has strongly anomalous velocity distribu- 
tions even in the absence of large-scale inhomogeneities. 
This is achieved by injecting the energy into rotational 
motion and allowing for a transfer to translations through 
strong friction. The system acts like a "transformer" con- 
verting Maxwellian degrees of freedom into distributions 
with fat tails. 
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